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Abstract 



A solid substrate, when exposed to a vapour, can interact with it in such 
a way that sufficiently close to liquid-vapour coexistence a macroscopically 
thick liquid wetting layer is formed on the substrate surface. If such a wetting 
transition occurs for a binary fluid mixture in the vicinity of the critical end 
point of demixing transitions, critical fluctuations lead to additional long- 
ranged interactions (Casimir forces) within the wetting layer, changing its 
equilibrium thickness. We demonstrate this effect by Monte-Carlo simulations 
of wetting layers of a symmetrical Lennard-Jones binary fluid mixture near its 
critical end point. The results suggest that the effect should also be detectable 
in corresponding wetting experiments. 
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I. INTRODUCTION 



Simple liquids or binary liquid mixtures in unbounded space are homogeneous and 
isotropic systems. In the vicinity of surfaces, however, these spatial symmetries are bro- 
ken and the fluid phase diagram exhibits new surface- induced features. Principal among 
these are surface phase transitions such as wetting and drying which can occur when a 
fluid comes into contact with a solid wall [Q . These transitions are brought about by long- 
ranged dispersion (Van der Waals) interactions between the particles of the fluid and those 
comprising the wall. 

Long ranged interactions of a quite different origin occur when a critical fluid is in contact 
with a wall |2|||. Criticality is characterised by a divergent correlation length and strong 
order parameter fluctuations. As a result of the long range correlations, the wall modifies the 
properties of the critical fluid over considerable distances and this can dramatically affect 
surface properties such as wetting behaviour. Thus, for example, a binary liquid mixture 
usually demixes close to a wall due to the preferential affinity of the wall material for one of 
the two components Q . If the demixing transition becomes critical, the adsorbed amount of 
the preferred component becomes macroscopic, a phenomenon known as critical adsorption 
II- 

If the system is constrained in more than one spatial direction, e.g. by the introduction 
of a second wall, the critical behaviour of the fluid is modified again. When the correlation 
length becomes comparable to the smallest linear extent of the system [0-0, the size de- 
pendence of thermodynamic quantities is governed by universal finite-size scaling functions. 
Such a finite geometry may be imposed eg. when fluids are adsorbed in slits and pores, or 
generated spontaneously in the course of a wetting transition eg. in a binary liquid mixture 
near its critical end point of demixing transitions |I|,[10|. 

An important consequence of confinement in near-critical fluids is the generation of long- 
ranged forces between the confining walls [11,12]. This phenomenon is a direct analogue of 
the well-known Casimir effect in electromagnetism [nj. Formally, the critical Casimir force 
arises as a result of the boundary conditions placed on the spectrum of order parameter 
fluctuations by the walls. This leads to a finite-size dependence for the singular part of the 



critical free energy, and hence a separation-dependent force between the walls [11,12]. Such 



critical Casimir forces are in addition to the usual dispersion forces which always exist in 
confined systems [P-l^|. However, they differ from dispersion forces in the major respect 



that they are governed by universal scaling functions |9j,ll,r^,^0[[. Precisely at the bulk 
critical point these scaling functions take universal values known as Casimir amplitudes 



11.2 



Much theoretical work has been directed in recent years towards an understanding of 
the Casimir effect in critical systems. For the Ising universality class, (which contains the 
critical demixing transition of binary liquid mixtures), exact results are available in d = 2 
dimensions for a strip geometry pT|j22"|| . In higher dimensions, exact results are limited 
to the Casimir force scaling function for the spherical model [O. For the experimentally 
relevant case of the d = 3 Ising universality class, only approximate results are available. 
These are based on real space renormalization |24| , the field theoretic renormalization group 
TT| , |20[] , and Monte-Carlo (MC) simulations p0|j25[] for the film geometry. Other geometries 



have also been investigated with regard to the possibility of performing direct surface force 
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measurements of the Casimir force using atomic force microscopy |26 , 27j (see also Ref. [|17|j). 

While the contemporary theoretical understanding of the Casimir effect in critical fluids 
is rather advanced, the experimental situation is not so satisfactory. Ideally, one should 
like to have experimental estimates of Casimir amplitudes and scaling functions to compare 
with theoretical estimates and to fill the gaps where no theoretical results are available. 
To date, however, no direct experimental demonstrations of the Casimir effect in critical 
fluid systems have been reported. This state of affairs seems to be attributable to the 
difficulties of performing high resolution surface force measurements within the extremely 
tight constraints of reduced temperature and pressure necessary to maintain criticality. In 
view of this, other more indirect routes to investigating the Casimir force have been sought. 
Perhaps the most promising of these is the effect of the Casimir force on wetting layers of 
a binary fluid mixture. If the vapour of a binary fluid is exposed to a solid substrate which 
strongly attracts the fluid particles, a thick liquid wetting layer will form on the substrate for 
state points sufficiently close to liquid-vapour coexistence. If, in addition, the temperature 
is chosen close to the critical end point (CEP) temperature (at which the line of critical 
demixing transitions intersects the liquid- vapour coexistence curve) , the presence of Casimir 
forces within the wetting layer will change its equilibrium thickness compared to a non- 
critical state point of the same undersaturation Since highly accurate measurements 
of wetting film thicknesses are possible using modern ellipsometry techniques, this would 
appear a rather promising approach to probing the Casimir force experimentally. As yet, 
however, there have been no clear reports of such an effect in experimental wetting studies 
of binary fluid mixtures. 

In view of the dearth of empirical results for the critical Casimir force, we have undertaken 
a computer simulation investigation of wetting behaviour close to the CEP of a binary fluid 
mixture. The aim of our study was two-fold: on the one hand to investigate the theoretical 
predictions of reference ||11|| , and on the other to try to make contact with real systems by 
working with as realistic a simulation model as feasible. Most previous simulation work on 
critical phenomena in confined geometry has taken the form of MC studies of lattice gas 
models, favoured for reasons of their computational tractability ^8,29|. Such simulations 



are usually performed using short-ranged surface fields to mimic the effect of walls. In this 
work however, we have employed an off-lattice Lennard- Jones (LJ) symmetrical binary fluid 
model at a planar solid wall, and have incorporated the effects of long-ranged dispersion 
interactions between the fluid and the wall. To our knowledge the wetting behaviour of such 
a system has only been previously studied for temperatures well below the consolute critical 
temperature |3(| . Hitherto, no attempts have been made to investigate directly the effect of 
criticality on the formation of wetting layers either in lattice or continuum fluid models. 

Our paper is organised as follows. Section || is devoted to providing some brief back- 
ground material concerning the bulk phase diagram of symmetrical binary mixtures and 
the theory of the effect of criticality on wetting layers. In section |111 A| we introduce our 
symmetrical binary fluid model and the MC simulation technique. We then detail our sim- 
ulation results for the wetting behaviour of the mixture, and the influence of the critical 
point. Finally, in section [TV] we summarise and discuss our results. 



II. BACKGROUND 
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A. Bulk phase behaviour of symmetrical binary mixtures 



The simplest model for a binary fluid mixture is a symmetrical mixture in which the 
pure components are identical, and only the interactions between unlike particle species 
differ. Such a model can exhibit a variety of phase diagram topologies depending on the 
relative strengths of the two types of interactions. In this work, we shall be concerned with 
symmetrical mixtures exhibiting a critical end point. The liquid- vapour phase diagram of 
such a system in the density-temperature (p — T) plane is shown schematically in figure |l|. 
At high liquid densities there exists a 'A' line of consolute critical points. Above this line 
the liquid is mixed, while below the line it is demixed. As the liquid density is decreased, 
the demixing transition temperature moves to lower temperatures. At some point, however, 
the liquid phase becomes unstable with respect to the vapour. The intersection between the 
A line and the liquid-gas coexistence curve marks the critical end point, which is unique as 
being the only point at which a critical liquid coexists with a non-critical gas. 

Critical end points have been the subject of some interest in recent years on account of 
the singularities they induce in the first-order coexistence phase boundary (in this case the 
liquid-gas boundary). These singularities have been studied both theoretically [[HJ and by 
simulation |32| . One is manifest as a bulge in the liquid branch density at the CEP, as shown 
in figure [l]. Another is found in the coexistence chemical potential. In the present work, 
however, we shall be concerned with the influence of the CEP on the wetting properties at 
the liquid-gas boundary and how it manifests the Casimir effect. In the following subsection 
we review existing theoretical predictions for the Casimir effect in wetting layers of a binary 
fluid mixture. 



B. Wetting layer thickness and criticality 

To analyse theoretically the wetting behaviour of the symmetrical binary fluid in the 
vicinity of the CEP, it is expedient to consider the phase diagram in the chemical potential- 
temperature (/i — T) plane, as shown schematically in figure ^|. Indicated on this diagram 
are the liquid- vapour coexistence line fi cx (T), the A line, and the CEP. 

Consider now the thermodynamic path labelled S in figure |2], which is given by the 
expression 

fJts(T)= f ji cx (T)-5 f jL, (2.1) 

where 5fi > is a constant. This path runs parallel to the liquid- vapour coexistence curve 
and passes the critical end point on the vapour side. In the presence of a sufficiently attractive 
wall (complete wetting regime), and for sufficiently small undersaturations (i.e. <5/i small), 
a thick liquid wetting layer will form at the wall for all points on S. An expression for the 
thickness of such a wetting layer is obtainable by considering the layer free energy per unit 
area (effective interface potential) as a function of its thickness /. In the limit of an infinite 
wall, this is given by (11]] 

= l(pi - p v )5fi + a wl + a tv + 8u(l) (2.2) 

Here the first term on the right hand side represents the free energy penalty of building 
up a liquid layer when the vapour is the stable bulk phase. The surface tension terms a w i 
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and o\ v contain the effects of short ranged interactions at the wall-liquid and liquid- vapour 
interfaces respectively, and are independent of /. The fourth term, Su(l), is a finite size term 
incorporating both the long ranged dispersion forces and any critical finite-size effects. 
Precisely at criticality, and for an infinite wall, 8u(l) takes the form [[11 



. /7N W k B TA . , 

-p+^p~- (2.3) 

In this expression, the first term on the right hand side represents the finite-size dependence 



of the (non retarded) dispersion forces [33 , whose amplitude is given by the Hamaker con- 



stant W. The second term represents the critical finite-size contribution to the free energy 



34j in three dimensions, controlled by a universal Casimir amplitude A, the magnitude and 
sign of which depend on the nature of the boundary conditions on the wetting layer. For A 
positive, the Casimir force is repulsive leading to a layer thickening, while for A negative, 



it is attractive resulting in thinning of the layer. The ~ in eq. |2.3| indicates that we neglect 
any density gradient in the liquid layer that may come about because of the attractive wall 
potential. A density gradient implies that not all the liquid layer is precisely at criticality. 
However, this is expected to result in relatively small corrections to eq. |2.3| . 

Minimising u(l) with respect to / yields the equilibrium layer thickness L c at criticality: 

V <wj - Pv) J 

For points on S away from the CEP, however, the critical finite-size term in eq. drops 
out. A useful quantity is thus the ratio of the equilibrium thickness at criticality L c to its 
value Lq away from criticality, for which one finds [|l2j 



L f A \ 1/3 

r i + k B T c £-) , (2.5) 



J C 



W 



to leading order in Sfi. Here we have ignored any temperature dependence of the Hamaker 
constant. 

This relationship expresses the change in thickness of a wetting layer as it becomes 
critical, and relates it to the Casimir amplitude A. As such it provides a potentially sensitive 
probe of the Casimir effect itself. 



III. MONTE CARLO STUDIES 



A. Model and simulation details 



The system we have studied is a symmetrical binary fluid, having interparticle interac- 
tions of the Lennard- Jones (LJ) form: 



u{rij) = Aeij 



12 



(3.1) 
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The following choice of parameters was made: a n = cr 22 = a 12 = = 1; e±x = £22 = e; 
612 = 0.7e. i.e. the pure components are identical, but the unlike interactions are weakened. 
In common with most previous simulations of LJ systems, the inter-particle potential was 
truncated at a distance of R c = 2.5cr. No long-range correction or potential shift was applied. 

The fluid was confined within a cuboidal simulation cell having dimensions P x x P y x D, 
in the x, y and z coordinate directions respectively, with P x = P y = P. As in previous work 
on LJ fluids |55], the simulation cell was divided into cubic sub-cells (of size the cutoff R c ) 



in order to aid identification of particle interactions. Thus P = pR c and D = dR c , with 
p and d both integers. To approximate a slit geometry, periodic boundary conditions were 
applied in the x and y directions, while hard walls were applied in the unique z direction at 
z = and z = D. The hard wall at z = was made attractive, using a potential designed 
to mimic the long-ranged dispersion forces between the wall and the fluid p6| : 



V(z) 













(t)1 



(3.2) 



Here z measures the perpendicular distance from the wall, e w is a 'well-depth' controlling 
the interaction strength, and we set a w — 1. No cutoff was employed and the wall potential 
was made to act equally on both particle species. 

Monte-Carlo simulations were performed using a Metropolis algorithm within the grand 
canonical (//, V,T) ensemble P7| , P3[ . Three types of Monte-Carlo moves were employed: 

1. Particle displacements 

2. Particle insertions and deletions 

3. Particle identity swaps: 1 —* 2 or 2 — > 1 



By virtue of the symmetry of the system, the chemical potentials [i\ and /x 2 of the two 
components were set equal at all times. Thus only one free parameter, // = \i\ = // 2 , couples 
to the density. The other variables used to explore the wetting phase diagram were the 
reduced well depth e/k B T and the reduced wall potential e w /k B T. During the simulations, 
the principal observable monitored was the total particle density profile: 

p(z) = [N 1 (z)+N 2 (z)]/(P 2 ) (3.3) 

which was accumulated in the form of a histogram. Here N\ and ./V2 are the number of 
particles of the respective species. Other observables monitored were the total interparticle 
energy and the wall interaction energy. 

The wetting layer thickness was obtained from the density profile as L = pj 1 J p{z)dz 
where p% is the liquid density. Although more sophisticated definitions of absolute layer 
thickness can be envisaged, this one is adequate for our purposes of detecting relative thick- 
ness changes arising from Casimir forces. 
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B. Method and results 



The bulk liquid-vapour coexistence curve properties of the symmetrical LJ fluid (with 
6x2 = 0.7e) have recently been investigated by MC simulation |32| . This work employed state- 
of-the art simulation techniques to investigate the nature of coexistence curve singularities 
induced by the critical end point. In the course of the study, highly accurate estimates were 
obtained for the location of the critical end point. Additionally, the locus of the liquid-gas 
coexistence curve fi cx {T) was measured to 5 significant figures in the neighbourhood of the 
CEP. 

Armed with accurate estimates for fi cx (T), one is in a position to perform detailed wetting 
studies very close to coexistence. Detailed knowledge of fi C x{T) is a prerequisite for obtaining 
a thick wetting layer and detecting changes in its thickness due to Casimir forces. Below we 
discuss the procedure we have employed for detecting such changes. 

Clearly for a liquid wetting layer to display critical behaviour at the bulk CEP, it must 
be sufficiently thick to exhibit quasi-3D properties. But for a thick wetting layer to form at 
a wall at all, the attractive wall potential must be sufficiently strong that complete wetting 
occurs at coexistence. To ensure that the CEP lay well within the complete wetting regime, 
a number of preliminary test runs were performed in which the temperature was set to 
its CEP value (T cep = 0.958) and the density profile p(z) monitored as coexistence was 
approached in a sequence of steps from the gas side. To achieve this, the chemical potential 
was simply increased towards its coexistence value in small increments of size A/x = 0.0025. 
The procedure was repeated for a number of different values of the wall potential e w . 

For small values of e„/^T < 1.7, it was found that the film thickness grew as coexis- 
tence was approached, but always remained finite (cf. figure 0(a)). The presence of a thin 
wetting layer right up to coexistence signals incomplete (partial) wetting. For stronger wall 
potentials (e^/Zc^T > 1.8), however, the film thickness became very large as coexistence 
was approached and ultimately almost completely filled the system (cf. figure |](b)). This 
suggests that the wetting transition lies in the range 1.7 < e w /ksT < 1.8. Nevertheless, 
to be quite sure that the CEP lay well above the wetting transition (ie. well within the 
complete wetting regime), all subsequent work employed e w /k B T = 3.0 at the CEP. 

We have studied the layer thickness L(fi S (T)) along the path /is(T) = fi cx (T) — 5fi. 
Given /z ca; (T), this path is specified solely by the value of which has to be chosen as 
a compromise between two competing factors. On the one hand, <5/i must be sufficiently 
small that a thick wetting layer forms. The thickness of this layer must be sufficiently large 
that quasi-3D critical behaviour occurs, and that the liquid-gas interfacial properties are not 
unduly influenced by short-ranged packing effect close to the attractive wall. On the other 
hand, one wishes to avoid having a layer that is very thick, lest the computational expense 
get out of hand. A minimum film thickness of L = 10a was judged sufficiently large for our 
purposes. However, since close to coexistence the layer thickness is extremely sensitive to 
changes in Sfi, fine tuning was necessary. It turned out that to obtain the required minimum 
thickness necessitated use of an extremely small <5/i, namely 5fi = 0.002. (corresponding to a 
relative undersaturation Sfi/p ~ 10~ 3 ). At this undersaturation, L was observed to fluctuate 
in the range 10cx ^ L < 13cr. 

Another condition to be satisfied in the simulations, is that the largest fluctuation in the 
layer thickness be small compared to the linear extent of the system in the z direction, ie. 
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L <<^ D. This latter condition ensures that the liquid gas interface does not interact with 
the hard- wall at z = D, which might otherwise obscure the effects of criticality, or even lead 
to capillary condensation. Accordingly the choice D = 40a was made. 

Ideally, in order to facilitate direct contact with real systems and existing theoretical 
work, one should like to simulate an effectively semi-infinite system having lateral dimensions 
P 3> L. In this limit, the finite-size critical behaviour of the system depends only on the 
ratio £/L of the correlation length to the layer thickness, and not on the lateral dimensions 
of the simulation box |3l |. Unfortunately, in this work we could not approach the semi- 



infinite limit as closely as we would have liked because of the huge computational expense 
entailed. In fact the largest lateral dimensions that we could reach were P = 12.5c, P = 15a 
and P = 17.5cr, containing average particle numbers of approximately iV = 1500, 1800, 3500 
respectively. We postpone discussion of the effects of the relative smallness of P until later. 

The limiting factor in the speed of the simulations and (hence the attainable system 
sizes) was found to be the extreme slowness of the interfacial fluctuations, which rendered 
it very difficult to accumulate accurate estimates of the average thickness. In an effort to 
ameliorate this problem, multi-histogram reweighting techniques were employed |38] . These 
allow one to combine data accumulated at individual state points and, by extrapolation, 
obtain estimates for observables at other not-too-distant state points. In our study the 
computational complexity meant that it was feasible to accumulate p(z) at just three points 
on the path S for each system size (see figure [| These simulation points were chosen to 
span the CEP temperature, namely T = 0.946, 0.958, 0.97. Subsequently, the data from the 
individual simulations were combined self-consistently and extrapolated to yield estimates 
for p{z) along the whole path. The data were collected from runs comprising in total 
approximately 3 x 10 9 Monte Carlo steps (MCS), where we define a MCS to be an attempt 
to perform each of the three types of MC move: particle displacement, insertion/deletion, 
and identity swap (cf. sec [III A|) . 

In figure f| we present our results for the layer thickness L(ns(T)), obtained by imple- 
menting this procedure. The main feature of the results for each system size P, is a peak 
in the film thickness close to the CEP. For the largest system size, the thickness either side 
of the peak is fairly constant, while for the smallest system size, the accessible range of 
temperature was not sufficient to encompass the whole peak. For the smallest system size, 
the critical thickening is ^ 10%, while for the largest system size it is ~ 5%. The peak 
position is at higher temperatures than the bulk CEP temperature and is closer to the CEP 
for the largest system size than for the smallest. The peak width also clearly narrows with 
increasing system size. We should caution, however, that the statistical quality of our data 
is not particularly high due to the difficulties mentioned above in collecting statistics. The 
smoothness of the temperature dependence for L(ns(T)) is to some extent an artifact of the 
histogram extrapolation procedure and we are not particularly confident of our estimates 
for the absolute film thickness. Nevertheless, we are confident of the ability of our procedure 
to identify changes in the wetting layer thickness. 

The results for the film thickening due to criticality can be compared with theoretical 
predictions for the semi-infinite system. In this limit, the change in layer thickness is given 
by eq. |2.5| . For our system, we calculate the Hamaker constant to be W ~ 2.5 at the CEP. 
The value of the Casimir amplitude A depends on the boundary conditions on the wetting 
layer. In the present model, these are of the form (+, 0) in the notation of reference [TTJ 
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This notation denotes an order parameter (concentration) that is pinned to a constant value 
at the wall due to the high particle density there, and which vanishes (on average) at the 
liquid-gas interface because of the low gas density. The most recent field theoretical and MC 
estimates |20| yield A(+,0) ~ 0.2. Inserting this into eq. |275| gives L c /Lq 1.025. Clearly 
this is a smaller thickening than we see, but the trend of our results is in the direction of 
this value and could lie close to the prediction for sufficiently large system size. 

In fact, it is possible to use finite-size scaling arguments to account for the direction of 
the observed trend in peak height with increasing P. To this end, let us rewrite the finite- 
size part of the critical effective interface potential |2.3| , taking account of the finite lateral 
extent of the system 



W 

MO = -p + 



k B T c A(+,0) 21A 



per 



I 2 



p3 



(3.4) 



Here we have simply added an additional scaling term to account for the finite system 
size P in the periodic x and y directions. The alteration to the free energy is assumed to 
scale like l/P, so that the contribution to u(l) (free energy per unit area) scales like l/P 3 . 
This new term can be interpreted as a finite-size correction to the chemical potential. 

Minimising this expression as before, gives the equilibrium critical thickness L c . The 
ratio of critical to non-critical thickness then follows as 



« 1 + 



k B T c A(+,0Y 
W 



1/3 



1 + 



2k B T c A 



per 



8fi(pi - p v )P 3 



-1/3 



The Casimir amplitude for periodic boundaries has been estimated as A 



per 



(3.5) 



-0.15 



25 



Thus the correction factor in P is larger than unity and decreases towards unity 
as P increases. If we take the liquid-gas density difference pi — p g ~ 0.5, one finds that 
the correction factor is approximately 1.06 for the P = 12.5a system size, reducing to 
approximately 1.02 for the P = 17.5cr system size. This accounts semi-quantitatively for 
the observed decrease in the size of the critical thickening as we increase P. 

The influence of the periodic boundaries is presumably also responsible for the the ob- 
served position of the peak in L(ps(T)). For a semi-infinite system the finite-size scaling 
function for ui{T) displays a peak for T < T cep and thus the peak in L(ps{T)) is also expected 
to occur for T < T cep . In our work, however, the peak in L(ps(T)) occurs for T > T cep . 
It seems reasonable, though, that were one to employ the finite-size scaling function for 
periodic boundary conditions in w(T) then this could actually change the temperature de- 
pendence of L(ps{T)) such that the maximum appears at some T > T cep , if P is not much 
greater than L. Unfortunately, the scaling function for a system with two periodic and one 
symmetry breaking boundary conditions is presently not known to any approximation. 



IV. CONCLUSIONS AND OUTLOOK 

In summary we have performed extensive Monte-Carlo simulations of wetting behaviour 
near the critical end point of a symmetrical binary fluid in contact with an attractive hard 
wall. The results demonstrate that critical fluctuations within the wetting layer engender 
Casimir forces, leading to small but significant changes in the layer thickness. As such, 
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they constitute the first empirical evidence for the Casimir effect in a critical system. Al- 
though an unambiguous determination of the limiting size of the critical thickening was 
hindered by large finite-size effects, the results are in order-of-magnitude agreement with 
theoretical predictions for a semi-infinite system. Indeed, a theoretical analysis of finite-size 
corrections to this scaling limit yields estimated size effects comparable to those observed 
in the simulations. Given the likely magnitude of the thickening in the semi-infinite limit 
(which we estimate to be ~ 3%), it seems likely that the effect should also be experimen- 
tally detectable in real (ie. non-symmetrical) absorbed binary fluids, many of which exhibit 
Hamaker constants of similar magnitude to the present model. 

With regard to the computational issues raised by our study, our results testify to the 
capability of modern simulation techniques to identify small changes in the wetting behaviour 
of realistic fluid models. Indeed, simulation is likely to prove of increasing value in the study 
of wetting by critical layers, since unlike field theoretical or density functional methods, it 
allows one both to tackle realistic systems and to deal properly with critical fluctuations. 
Nevertheless, algorithmic improvements are clearly desirable and necessary if one is to realise 
geometries that are effectively semi-infinite. The chief problem experienced in this work 
was of extremely slow fluctuations of the liquid-gas interface. In view of this, it may well 
pay dividends in future work to employ a simulation algorithm that focuses more of the 
computational effort on the interfacial region itself, this being the bottleneck for phase 
space evolution. Such algorithmic improvements are the matter of ongoing work. 

Finally, looking ahead to future work, it would doubtless be interesting to investigate the 
nature of the interfacial properties between the critical liquid and non-critical gas phases. 
This matter has already been the subject of a number of detailed theoretical investigations 
aimed at elucidating the temperature dependence of the interfacial shape and surface tension 
on the approach to the CEP P^fTO||?T| |. Given present capabilities, a simulation study of 
interfacial properties at the CEP would certainly seem feasible and would nicely complement 
existing experimental studies [O, 
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FIGURES 




FIG. 1. Schematic representation of the phase diagram of a symmetrical binary fluid mixture in 
the density-temperature plane. The full curve is the liquid-gas coexistence envelope. The dashed 
curve is the A-line of critical demixing transitions. The two curves intersect at the critical end 
point. The singularity in the liquid branch at the CEP is also shown p2| . 
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Temperature T 



FIG. 2. Schematic phase diagram of a symmetrical binary fluid. The full curve is the first 
order liquid-gas phase coexistence line [i cx (T). The dashed curve is the critical (A) line of second 
order transitions separating the mixed and demixed liquid phases. The two curves intersect at the 
critical end point. Also shown is a path S parallel to fi cx (T) on the gas side. 
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FIG. 3. (a) Density profiles p(z) on the approach to coexistence for e w = 1.7. The film 
thickness remains small right up to coexistence signifying incomplete wetting, (b) Density profiles 
for e w = 1.8. The film thickness grows very large as coexistence is approached, signifying complete 
wetting. 
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FIG. 4. The thickness of the wetting layer as a function of temperature along the thermody- 
namic path <S defined in the text. Data is shown for each of the three system sizes studied. The 
results were obtained from multihistogram extrapolation of simulation data accumulated at three 
points on this line, corresponding to temperatures T = 0.946, 0.958, 0.97. 
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